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Using a two dimensional lattice model we investigate the crack growth under the influence 
of remote tensile forces as well as due to an internally applied pressure (hydraulic fracturing). 
For homogeneous elastic properties we present numerical finite-size scalings for the breaking 
stresses and pressures in terms of crack lengths and lattice sizes. Continuum theory predicts 
for the tensile and for the pressure problem identical scaling functions. Our findings for the 
tensile problem are in very good agreement with continuum results. However, for the hydraulic 
fracture problem we observe a different finite-size scaling. We explicitly demonstrate that the 
modified scaling is a consequence of the discrete structure of the lattice (micro structure). 
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I. INTRODUCTION 



About 30 years ago considerable interest has been 
devoted to the properties of generalized Cosserat elas- 
ticity. One of the main activities was the formulation 
of continuum theories reproducing lattice dispersion 
relationships to a certain extent [jlj . These approaches 
(gradient theories) were designed in order to incorpo- 
rate certain atomistic aspects into generalized elastic 
continua. However, it was early recognized that such 
theories were in principle also applicable to the de- 
scription of micro structural and dislocation phenom- 
ena Q . Recently a first order (strain) gradient theory 
was applied to a fracture mechanical problem obtain- 
ing results reminiscent of Barrenblatt's cohesion the- 
ory fill- However, the gradient elasticity approach 
leads to, in general, formidable mathematical diffi- 
culties, to a somewhat suspect choice of additional 
boundary conditions, and to numerous physically not 
well understood new material constants. 

These difficulties are circumvented in discrete frac- 
ture models which have been adopted during the past 
decade in order to study the effect of structural het- 
erogeneities and fluctuations on fracture growth pro- 
cesses under various physical conditions. Using tri- 
angular networks of springs stretched on the outer 
boundary into the six directions of a hexagon the 
breaking of a material from a central hole has been 
investigated by several authors pj. Employing disor- 
dered beam networks, the scaling of micro crack den- 
sities Q and the growth of cracks due to an inner 
hydraulic pressure Q have been studied in some de- 
tail. These investigations offer a new methodological 
approach to study failure processes in heterogeneous 
brittle materials on a mesoscopic length scale close 
to experimental situations [|| . Our interest is focused 
on brittle materials because their rheological behavior 
can be regarded most simple as they do not show any 
remanent (plastic) deformations ||. 

One difficulty arising during the examination of nu- 
merically observed results often consists in finding ap- 
propriate scaling variables. As we will see below this is 



not a trivial task, not even for perfect straight cracks. 
We focus on finite-size scaling because the employed 
lattice sizes L are always finite one is, however, mostly 
interested in properties of very large or infinite sys- 
tems. Thus the finite-size scaling approach appears 
to be the only conclusive numerical method to verify 
(conjectured) asymptotic behavior for L — > oo. 

In the present paper we will show that the micro 
structure (lattice structure) can change both asymp- 
totic and finite-size scaling of the breaking characteris- 
tics, depending on the employed boundary conditions. 

In Sec. H we briefly outline the fracture model, the 
employed boundary conditions and how to determine 
the external breaking str esse s and pressures. Follow- 
ing this we present in Sec. HI the main finite-size scal- 
ing results for breaking stresses and pressures of linear 
cracks on two-dimensional beam lattices. We compare 
our findings with known analytical continuum results. 

For crack growth scaling properties in disordered 
materials the reader is referred to Ref. |(| 



II. THE MODEL 

We investigate both the case of uniaxial loading 
due to remote tensile forces and the loading due to 
a hydrostatic pressure acting from the interior on the 
crack surface. The former situation is typically en- 
countered when a material is pulled at its outer faces, 
while the latter one can be found in engineering appli- 
cations of hydraulic fracturing We consider the 
beam model on a two dimensional square lattice of 
linear size L J8|,[ll| . Each of the lattice sites i carries 
three real variables: the two translational displace- 
ments Xi and yi and a rotational angle tpi. Neigh- 
boring sites are rigidly connected by elastic beams of 
length I. The beams all have the same cross section 
and the same elastic behavior, governed by three ma- 
terial dependent constants a = l/(EA), b = l/(GA) 
and c = I 3 /(EI) where E and G are the Young 
and shear moduli, A the cross section of the beam 
and / the moment of inertia for flexion. We used 



for all simulations the values a = 1.0, b — 0.0017 
and c = 8.6. When a site is rotated (ipi ^ 0) the 
beams bend accordingly always forming tangentially 
right angles with each other. In this way local mo- 
menta are taken into account. For a horizontal beam 
between sites i and j one has the longitudinal force 
acting at site j: Fij = a(xi — Xj); the shear force: 
Sij = l3(yi — Uj) + ^KVi + fj)i an d the flexural torque 
at site j: M y - = f - yj + lipj) + Sl 2 {(pi - (fj), using 
a = 1/a, (3 = 1/(6 + c/12) and S = /3(6/c+l/3). Cor- 
responding equations hold for the vertical beams. The 
hydrostatic stress acting on a beam is given by its lon- 
gitudinal force per beam cross section, <7y = Fij /A. 
In mechanical equilibrium the sum over all internal 
and external forces (torques) acting on site j must 
vanish. We will neglect inertial forces, assuming a 
sufficiently slow fracture process. The cohesive prop- 
erties of the beams are considered to be homogeneous, 
i.e. the beams behave linear-elastic up to their the- 
oretical strength F t h (cohesion force), above which 
they irreversibly break (zero elastic constants). At 
the beginning of each simulation we break one ver- 
tical beam located at the center of the lattice. This 
is the initial crack. Corresponding to the employed 
boundary conditions described below we calculate the 
internal force distribution Fij of beams connecting 
sites i and j using a conjugate gradient method ]l^ |. 
From this longitudinal forces we determine and break 
the most over stressed beam, i.e. the beam carry- 
ing the highest force, Fi j = maxjjjj. F^. The max- 
imum force allows us to calculate a scaling factor 
A = F t h/ Fi j = o'th/o'iojo which in turn is used to de- 
termine the macroscopic breaking stress a c or break- 
ing pressure P c (measured in units of oth) necessary 



to fulfill the local fracture criterion er.; 



a th . Such 



scaling is possible due to the linearity of the employed 
elastic equations. Breaking the beam connecting sites 
io and jo destroys the balance of forces at those sites 
and one has to relax the system to its new equilibrium 
configuration. Then the above steps are repeated until 
the lattice breaks apart. 

It is well known that the highest stress enhance- 
ment factors occur at the crack tips for the pres- 
sure as well as for the tensile problem. Because we 
consider only homogeneous cohesion forces the terms 
'most over stressed' and 'highest stress' are equivalent 
and the cracks grow linear in direction of the highest 
local tensile force. In our case the direction of crack 
growth is the x-direction because of the uniaxiality of 
imposed boundary conditions, see below. 



A. Tensile Fracturing 

To impose an external strain we attach at the bot- 
tom of the lattice a zeroth line on which for all sites 
x i — Hi = ¥>i = are fixed, and on the top we attach 
a (L + l)st line on which all sites have the same fixed 
values Xi = ipi = and yi = 1. With this boundary 
conditions the external displacement in y-direction is 
fixed to unity (Dirichlet boundary conditions) and one 
can imagine them as being represented by rigid bars 
attached at the top and the bottom of the lattice . 



For the boundary conditions in x-direction we con- 
sider periodic as well as free lattice boundaries. The 
externally applied stress, Oo, necessary to maintain a 
unity displacement between the two rigid bars can be 
easily calculated by summing up all yi displacements 
over the first line, Co = Uii where a is the lon- 

gitudinal force constant, A the cross section of the 
beam and L the (dimcnsionless) lattice size. Employ- 
ing the above defined scaling factor A we express the 
externally applied breaking stress, cr c , in the form, 



= Actq, 



or in dimcnsionless form, 



o_c_ 



1 a ^ 



F, 



(2.1) 



(2.2) 



With this notation the external breaking stress a c is 
just the stress necessary to yield the local breaking 
condition <7;„ 



Oth- 



We will use Eq.(2.2) for our finite-size scaling pur- 
poses. It expresses the breaking stress, measured in 
units of cohesion stress, in terms of the beam force 
Fi j at the 'crack tip' and the average force ^ ^2 Hi 
on a beam due to the global unit displacement. In 
general both terms depend on the number of broken 
beams N (crack size) and on the lattice size L, see 



Sec. HI A 



B. Hydraulic Fracturing 

Experimentally hydraulic fracturing is encountered 
when an incompressible fluid is injected under high 
pressure into an existing crack. The characteristics of 
this particular loading mode is that the loading of the 
crack only happens on the crack surface itself, but not 
on remote surfaces like in the case of tensile fractur- 
ing. Recently hydraulic fracturing has been numeri- 
cally addressed by introducing for each broken beam 
a force dipole of strength Fq simulating an inner pres- 
sure Pq = Fq/A [0. In the present paper we will 
follow this implementation. Starting from one verti- 
cal broken beam (zero elastic constants, one vertical 
force dipole) we calculate the internal stress distri- 
bution F^ using the conjugate gradient. Calculating 
the scaling factor A as described above we determine 
and break the beam carrying the highest stress, Fi j . 
Thereafter the lattice has to be relaxed again (now two 
broken beams, two force dipoles) and one repeats the 
above mentioned steps until the lattice breaks apart. 

We have performed calculations for free external 
boundaries in x- and y-direction as well as for peri- 
odic boundaries in x-direction and free boundaries in 
y-direction. As we only impose stresses the elastic 
solution (displacement field) is unique up to a gen- 
eral translation and rotation. Without restriction we 
therefore fix the displacement of an arbi trar y lattice 
site to zero. Analogous to Eqs. ( |2.1[ )and (2.2) one has 
for the breaking pressure, 



Pr = APn 



(2.3) 



and 



3l 

o~th 



40J0 



(2.4) 



respectively. It should be noted that o~th is the tensile 
strength and not the compressional strength. In fact 
the stresses at the tips are tensile stresses as the crack 
is hydraulically ope ned. 

We will use Eq. (2.4) for the finite-size scaling of 
the hydraulic problem. 



III. RESULTS 

As explained in Sec. || we obtain from our simula- 
tions the macroscopic breaking stresses and pressures, 
<j c {N, L) and P C (N, L), measured in units of the theo- 
retical strength a t h , m terms of the number of broken 
beams N and the employed lattice sizes L. As we are 
interested in the finite-size scaling properties of the 
breaking stresses (E c = a c ) and pressures (S c = P c ) 
we introduce two new functions, iff(N) and $>(N/L), 



V c {N,L)=<j th y(N)<S>(N/L). 



(3.1) 



The function $(JV/L) describes the finite-size cor- 
rections to the scaling iff(N) between the breaking 
stress/pressure and the crack length for given bound- 
ary conditions in the infinite system. The quantity 
N/ L is for our purposes the most suitable scaling vari- 
able as it represents the length of a linear crack mea- 
sured in units of the employed lattice size. In general 
both functions, iff and $ will depend on the employed 
boundary conditions. In order to show numerically 
that the finite-size scaling function Q(N/L) exists, one 
needs to know (or to conjecture) the explicit form of 
iff(N), i.e. the asymptotic behavior of S c as L — > oo. 

We will access the latter information from sim- 
ple continuum mechanical considerations of symmet- 
ric elasticity. 



A. Tensile Fracturing 

Previous continuum mechanical scaling considera- 
tions have successfully been applied to fracture mod- 
els. It has been shown jl3| that the breaking stress 
o~ c for tensile fracturing of a finite, periodic beam lat- 
tice follows very well the known 'Tangent-Formula' of 
fracture mechanics Jl4|,[l5| , 



CTth 



N' 



-1/2 



N ttN 

r ot u' 



(3.2) 



By comparison with Eq. (3T) we see that the asymp- 
totic scaling for the breaking stress a c is given by 
v]/ f« N^ 1 / 2 , a result that originally was found by Grif- 
fith. In Fig. we show the finite-size scalings for the 
numerically obtained tensile breaking stresses (a) for 
free boundary conditions and (b) for periodic bound- 
ary conditions. For both curves the data collapse is 
quite acceptable. While for the periodic pr oble m the 
finite-size correction $(N/L) follows Eq. ([T|]) Q, 
no analytical expression for <E> is known in the case 
of free boundary conditions. However, the differences 
for the breaking stresses at given values of N/L are 



less than 0.1<j t h, see Fig. [|, and Eq. (3.2) represents 
a reasonably good approximation for finite solids for 
many practical purposes. 



B. Hydraulic Fracturing 

In general analytical solutions of crack problems 
are difficult to obtain and even for two-dimensional 
problems in infinite, semi-infinite or periodic domains 
one is often faced with complicated integral equations. 
The reader interested in this is referred to Muskhel- 
ishvili @ and Sneddon 

It has been shown analytically that the breaking 
pressure P c of a semi-infinite two-dimensional contin- 
uum with periodic boundary conditions is given by 
Eq. ( ft.2| ) where one has to replace a c by P c 
show in Fig. || the analogous plot to Fig. 
breaking pressure. One does not obtain a data col- 
lapse, indicating that iff ~ iV -1 / 2 , see Eqs. (3.1) and 
( |3.2| ) , is not the scaling of the breaking pressure in the 
asymptotic limit of infinite lattices. 



1§. We 

for the 



This result is somewhat unexpected, because the 
continuum limit predicts identical scaling relations for 
breaking stresses and pressures and on the other hand 
we find good agreement between lattice and contin- 
uum scalings for the tensile problem. 

However, the proposed model for hydraulic fractur- 
ing is defined on a lattice and not in a continuum. As 
lattices always show additional structure, higher order 
gradient terms of the continuum displacement field 
appear in the elastic solution, which have no counter- 
parts in simple continua |jl7|| . 

In the following we present an argument why the 
lattice finite-size scalings for the pressure and tensile 
problem are different. 

In the tensile problem the loading forces are act- 
ing at the lattice boundaries remote to the stress free 
crack surfaces. Therefore local stresses close to the 
crack tips are carried through nearly the whole elas- 
tic volume. Contrary to this in the hydraulic problem 
loading only happens due to force dipoles acting at the 
crack surfaces close to the crack tips. One therefore 
might expect that the breaking pressures P c show a 
more pronounced '(micro structure) lattice behavior' 
than the breaking stresses a c . For very large cracks 
these differences should diminish as the discrete load- 
ing approaches a force density. 

In order to extent a linear crack of length N in 
an infinite continuum one has to impose a breaking 
pressure P c ~ iV -1 / 2 . However, employing a single 
double force of magnitude F c one finds, F c ~ N 1 / 2 for 
the critical loading force (related Boussinesq problem) 
©■ 

We argue now that for a single crack in an infinite 
disk loaded by equidistant force dipoles the asymp- 
totic scaling for iff is equivalent to that of an infinite 
lattice. Knowing iff we plot the finite-size scalings for 
the investigated lattices. 

The asymptotic continuum scaling for the breaking 
pressure (as calculated from equidistant force dipoles) 
is most easily obtained using the complex stress func- 
tion of Westergaard [Q^2|. We find for the asymp- 
totic scaling in an infinite disk, 



*(JV) 



ATl/2 



2 1 + Nf(N)'' 



with 



(N-l)/2 



fc=i 



(jV±i)2_ fe 2 



(3.3) 



AT odd. (3.4) 



The sum f(N) converges rapidly towards arcsin(-^— j-) 
and for large N towards ir/2. Hence in the limit of 
a large number of force dipoles we obtain the afore- 
mentioned asymptotic Griffith scaling for the breaking 
pressure, P c /<Jth = iV~ 1/2 . 

In Fig. U we show the resulting finite-siz e sca lings 
for the breaking pressure, based on Eqs. (3.1) and 
( |3.3| ). In this plot a reasonable scaling is obtained, 
which proves our earlier statement that the asymp- 
totic limits (L — > oo ) of the breaking pressures for 
continua and lattices are different. Comparing Fig. ^| 
with Fig. [j] one also notes a significant deviation from 
the scaling form Eq. ( |3.2| ), demonstrating a modified 
correction to scaling behavior. 



IV. CONCLUSIONS 

We have presented finite-size scaling results for 
breaking stresses (tensile fracturing) and breaking 
pressures (hydraulic fracturing) of elastically and co- 
hesively homogeneous beam lattices in two dimensions 
(plane stress condition) . For the tensile problem it has 
been shown that both the asymptotic lattice scaling 
^(N) as well as the finite-size lattice scaling §{N / L) 
are well represented by the known continuum scaling 
form ('Tangent-Formula'). 

In contrast the hydraulic fracture problem formu- 
lated on a lattice, has lead us to different asymp- 
totic and finite-size lattice scalings, though the cor- 
responding continuum scalings for breaking pressures 
and stresses are identical. Explicitly we constructed 
the asymptotic lattice scaling of breaking pressures 
from purely continuum mechanical considerations, by 
taking into account the discrete nature of the loading. 

As the overwhelming number of brittle materials is 
granular in character, the interesting question arises 
whether their breaking characteristics can be satis- 
factory described by continuum models regardless of 
employed loading conditions. 
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FIG. 1. Finite-size scaling for the rescaled breaking stress a c / utuN 1 ^ 2 as a function of the rescaled crack length N/L, 
(a) for free boundaries in x-direction and (b) for periodic boundaries in x-direction. Lattice sizes: (A) for L = 20, (+) 
for L = 40 and (o) for L = 100. 
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FIG. 2. Plot of the rescaled breaking pressure Pc/crthN 1 ^ 2 as a function of the rescaled crack length N/L for periodic 
boundaries in x-direction. In contrast to Fig. [I] we do not find any data collapse, indicating that ^(N) ~ TV -1 / 2 in 
Eq. ( |3.l[ ) is not the appropriate asymptotic expression for the breaking pressure. Lattice sizes: (A) for L = 20, (+) for 
L = 40 and (o) for L = 100. 
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3. Finite-size scaling for the rescaled breaking pressure P c /ath $ 1 {N) with <I> 1 (jV) = ^ 1+ ^//2 N ^ ■. 



as a function of the rescaled crack length N/L; (a) for free boundaries (A) : L = 50, (+) : L = 100, (o) 
for periodic boundaries (A) : L = 20, (+) : L — 40, (o) : L = 100. Note the data collapses. 



see Eq.flyD, 
L = 300; (b) 



